Spatially dependent Rabi oscillations: 
an approach to sub-diffraction-limited CARS microscopy 
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We present a theoretical investigation of coherent anti-Stokes Raman scattering (CARS) that is 
modulated by periodically depleting the ground state population through Rabi oscillations driven 
by an additional control laser. We find that such a process generates optical sidebands in the CARS 
spectrum and that the frequency of the sidebands depends on the intensity of the control laser light 
field. We show that analyzing the sideband frequency upon scanning the beams across the sample 
allows one to spatially resolve emitter positions where a spatial resolution of 65 nm, which is well 
below the diffraction- limit, can be obtained. 

PACS numbers: 31.15.ap, 87.64.-t, 78.20.Bh, 78.47.Fg 



I. INTRODUCTION 

Optical microscopy is one of the key techniques used 
to analyze biological processes at systemic, cellular, and 
sub-cellular levels. However, current optical microscopy 
techniques have insufficient resolution to observe interac- 
tions at the resolution desired by cell biologists-that is 
at the molecular or single functional group level. 

Current near-field techniques are sub-diffraction- 
limited, such as scanning near-field optical microscopy 
[l[ but they suffer from the disadvantage that they are 
limited to surfaces. Two notable far-field techniques 
have achieved sub-diffraction-limited resolution: stim- 
ulated emission depletion (STED) microscopy 0], and 
stochastic optical reconstruction microscopy (STORM) 
0, 0] ■ Both techniques require that the sample is labeled 
with fluorescent dyes. In addition, STED requires that 
the intensity at the sample is rather high to resonantly 
saturate the electronic transitions of the label, hasten- 
ing photo-bleaching and potential cyto-toxicity 0, Q. 
STORM, on the other hand, uses very low light inten- 
sities, but requires very special labels that can be inten- 
tionally switched between dark and light states. Further- 
more, STORM images take a substantial amount of time 
to acquire obstructing live cell imaging. 

Ideally, biologists would prefer to employ label-free 
imaging techniques, such as CARS microscopy, which re- 
places labeling by detecting the chemically specific vi- 
brational modes of the molecules. Furthermore, the non- 



resonant nature of the excitation process limits photo- 
damage [2j . Techniques to achieve sub-diffraction-limitcd 
CARS imaging are the subject of active investigation. 
Linear techniques to suppress CARS emission have been 
the subject of both experimental and theoretical inves- 
tigation Q. However, linear techniques cannot be used 
to obtain sub-diffraction-limited resolution. In previous 
work, we investigated a saturation process, analogous to 
STED, that prevents the build-up of the vibrational co- 
herence required for CARS emission [9]. However, the 
saturation was presumed to be fully incoherent, via a 
combination of long lifetimes and short dephasing times, 
leading us to explore alternative approaches that also in- 
clude coherent effects. 

Here, we use a density matrix approach 10] to model 
the CARS emission process. We consider CARS emis- 
sion in a four- level system (see Fig. [IJ in combination 
with an additional light field that couples an additional 
level |4), called control level, to the ground-state level |1). 
A mechanism that leads to sub-diffraction-limited CARS 
microscopy has been identified for this system: for the 
case that the lifetime of |4) and the dephasing of the 
the |1)-|4) transition are long compared to the inverse 
transition rate imposed by the additional laser, Rabi 
oscillations between |1) and |4) induce a Rabi-splitting 
of the CARS emission. The Rabi-splitting is intensity- 
dependent, which can be used to resolve features within 
a diffraction-limited volume. 



II. THEORETICAL FRAMEWORK 
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A detailed description of the level scheme, density 
matrix equations, and light fields has been given else- 
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FIG. 1: Energy level diagram for the CARS process with an 
additional level |4). Level |1) is the ground level and ini- 
tially fully occupied. Level \2) is a vibrational level of the 
medium, |3) is the excited level and level |4) the control state, 
which has a low decoherence rate. The arrows indicate possi- 
ble transitions induced by the incident laser fields, which are 
far detuned from the transition frequency (uj p , u>s and ui pr ), 
or on resonance (lo c ). Through coherent anti- Stokes Raman 
scattering (CARS), the medium emits light at the frequency, 
wcars, while the additional control field, uj c , drives Rabi os- 
cillations at the |1) - |4) transition. 



where Q • In summary, the four-level system (see Fig. [T]) 
has a ground state (|1)), a vibrational state (|2)), an ex- 
cited electronic state (|3)), and the control state (|4)). 
Transitions between the ground and excited states, the 
vibrational and excited states as well as the ground and 
control states are dipole allowed, while all other tran- 
sitions are dipole forbidden. The medium is irradiated 
with four pulsed light fields, two of which (pump, oj p ; 
Stokes, cos) are two-photon resonant with the |1) - |3) 
transition, driving population into |2) and inducing a 
coherence between |1) and |2) (we refer to this as the 
vibrational coherence). The vibrational coherence in- 
duces two sidebands on the probe light field (to pr ), one 
of which is known as the CARS emission frequency 

(iOp — ids + w pr — ^cars)- 

The additional laser (called the control laser uj c ), res- 
onant with the |1) - |4) transition, influences the CARS 
emission. To understand the control laser's effect, we cal- 
culate the temporal evolution of the laser pulse envelopes 
and the envelopes of the density matrix elements, along 
with their radiating fields, as described in detail in ref- 
erence We study CARS emission for the case where 
|4) has a decoherence rate that is slow compared to the 
duration of the laser pulses, resulting in Rabi oscillations 
between |1) and |4). Although |4) is considered to be a 
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vibrational state in this calculation such as in 
any dipole allowed state with the requisite decoherence 
rate can be used. 

For generality, we choose energy levels, detunings and 
pulse durations that are typical for CARS emission pro- 
cesses from molecules. The |1) - |3) transition angu- 
lar frequency is set to 4,700 THz (approx. 400 nm or 
25,000 cm" 1 ), and the |1) - |2) frequency to 314 THz 
(approx. 6 /zm or 1,667 cm -1 ). Likewise, we choose a 
|1) - |4) transition frequency of 565 THz (approx. 3.33 /im 
or 3,000 cm -1 ). 

The combination of the pump laser wavelength as 



800 nm, and Stokes as 923 nm is taken to provide two- 
photon resonance with the |1) - |2) transition. The probe 
laser wavelength is taken as 600 nm. The total lifetime 
of state 1 3) and the lifetime of state |4) are taken to be on 
the order of nanoseconds [13| , while the decoherence rates 
between states involved in the CARS emissionprocess are 
of the order of picoseconds (see e.g., [H, H, HB| ) • All the 
laser pulse durations r are set to 7 ps (1/e 2 ) except for 
the control laser field, which is continuous during the in- 
terval of CARS emission. The numerical calculations ex- 
tend over 30 ps in steps of 0.1 femtoseconds. The amount 
of output data was set to 70,000 points for fast compu- 
tation and causes only a small amount of discretization 
noise. 



III. RESULTS 

We performed a series of calculations of the emission 
spectrum around the CARS frequency. In these calcu- 
lations lu c was taken to be resonant with the |1) - |4) 
transition (the detuning A14 = 0), while the decoher- 
ence rate of this transition was set to Ti4 = 0.1 THz. 
The amplitude of the control field and, hence, the Rabi 
frequency (Qr) were varied in steps. The results are 
shown in Fig. [2] where the spectrum around the CARS 
frequency, as emitted by the medium, for a number of 
values of Qr is given. A single emission peak is seen at 
low Rabi frequencies and indicates standard CARS emis- 
sion, as expected. It can be seen that, for increasing Qr, 
the spectrum shows two symmetrically placed sidebands. 
From our calculations, we find that these sidebands show 
a spacing from the CARS carrier frequency that coincides 
with Qr, given by [l(|: 
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Here E c = \J (2/ c /eo) is the electric field of the control 
laser and I c the control intensity, and fiu is the dipole 
moment at the control transition. From this, we con- 
clude that the Rabi oscillation periodically depletes the 
ground state, which, in turn, modulates the vibrational 
coherence. As a result, the CARS emission becomes am- 
plitude modulated which spectrally shows up as two Rabi 
sidebands. The splitting of several THz is large enough 
to be distinguishable in a real CARS experiment using 
standard spectrum analyzers. 

It is critical, for obtaining a noticeable Rabi splitting in 
the CARS spectrum, that the frequency of the Rabi oscil- 
lations exceeds the decoherence rate Ti4, as can be seen 
in Fig.rj^b), where the spectrum of the CARS emission, 
modulated with Qr = 4 THz, is shown for increasing 
values of the decoherence rate. Up to decoherence rates 
which correspond to sub-picosecond decoherence times, 
the 4 THz modulation is still detectable in Fig.^b). Ob- 
taining such high Rabi frequencies is important, because 
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FIG. 2: (a) Emission spectra containing the CARS emission 
peak for increasing values, on a logarithmic scale, of the in- 
duced Rabi frequency oscillations between states |1) and |4). 
The amount of splitting of the sidebands is equal to Qr. (b) 
Emission spectra containing the CARS emission peak for in- 
creasing values, on a logarithmic scale, of the decoherence 
rate of Fi4 with Qr — 4 THz for each spectrum. Note 
that the amplitude of the sidebands decreases and disappear 
around where Qr equals the decoherence rate. The increase 
in Rayleigh scattering of the incident light fields for the lower 
set of decoherence rates, leads to an increase in background 
signal at the CARS emission wavelength. 



the decoherence times in liquids are very short, typically 
on the order of 5 ps [ll[ • 

The induced sidebands can be used to obtain positional 
information by noting from Eq. [T]that the Rabi splitting 
of the sidebands depends on the intensity of the con- 
trol laser. Furthermore, the geometrical alignment of the 
dipole moment /ii4 may be obtained by rotation of the 
polarization angle of the control laser beam, providing 
additional information on the structure of the sample. 
Considering a single emitter illuminated by a Gaussian 
beam, we note that the radial distance of the emitter from 
the center of the control beam can be obtained by mea- 
suring the exact frequency of the sidebands, which can 
be called a Rabi-labeling of the emitter position. Such 
a sp ace-dependent Rabi frequency is a well-known effect 
[161 ] . however, it is usually considered to cause undesired 
broadening and is eliminated, e.g., in pump-probe exper- 
iments, by using a much larger pump beam as compared 
to the probe beam 17 1. 



The absolute position of an emitter location can be 
calculated by trilateration. We illustrate this with a two 
dimensional example, however, a similar approach also 
applies to three dimensions. Consider an emitter located 
at the Cartesian coordinates x\ and y±. The observed 
CARS emission spectrum, with its sideband frequency 
depending on the local intensity of the control laser, will 
not determine these coordinates but instead provides the 
distance from the control beam center, r\ = y 'x\ + yf, 
which is a ring of possible locations centered around the 
control beam. By scanning the control laser a known 
distance dx along the x axis, the radial location of the 
emitter is changed to x\ - dx, y\. As a result, the emit- 
ter modulation frequency will c hange, revealing a new 
distance, r2 = y (xi — dx) 2 + y\. The expressions for r\ 
and ri can then be solved for xx, and a subsequent scan 
along y by a known dy determines y\ as well. Turning 
our attention to the spatial resolution of such a position- 
ing scheme, we show that the scheme can provide CARS 
with sub-diffraction-limited resolution. 

We note that the control laser, resonant with the 
1) - |4) transition and focused to a diffraction-limited 
Gaussian-shaped intensity profile, centered on the other 
laser beams, illuminating a distribution of emitters, will 
generate a distribution of Rabi sidebands on the CARS 
signal. If the control laser has a peak envelope field 
strength of Eq and a 1/e-spot-size of wq, then the range of 
observed Rabi frequencies within the diffraction-limited 
spot is Umax = \h*iaEq / h\ m the center to £l m in = 
1 / e \n14E0 / h\ at radial distance wq. The radial reso- 
lution is then determined by how accurately the Rabi 
splitting frequency can be measured in the CARS spec- 
tra. Given a Rabi frequency, f2, corresponding to a radial 
distance from the center of the focus of r, and a frequency 
measurement accuracy of d£l, corresponding to a radial 
resolution of dr, it can be shown that 
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It can be seen from the transcendental equation [2] that 
dr becomes smaller than the diffraction-limit (smaller 
than wq), and that higher Rabi frequencies lead to higher 
resolution, while emitters located away from the center of 
the control beam are resolved better than emitters close 
to the center. The resolution is expected to be maxi- 
mal at the steepest intensity slope of the control beam, 
therefore, by using a Gaussian shaped control beam, the 
highest resolution is expected to be found off center. This 
is illustrated in Fig. [3] where we have plotted the Rabi- 
labeling resolution as a function of both Rabi frequency 
and radial location, for dfl = 3 cm -1 (which is a typi- 
cal value for the linewidth of the CARS peak in a spec- 
trum). This figure shows that considerable improvement 
over diffraction-limited CARS microscopy is possible. 

To make the improvement explicit, we consider a spe- 
cific example for a microscope with a numerical aperture 
of 1.2 and laser wavelengths that were used in our cal- 
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FIG. 3: Spatial resolution as a function of the Rabi frequency 
and radial location of the emitter. 

culations. In this case, the diffraction-limited resolution 
is 171 nm. If the control laser (A=3.3 /jm, focused to a 
diffraction-limited spot, induces 100 cm -1 Rabi oscilla- 
tions, then at a radial distance of 1.2 /xm, the resolution is 
65 nm. We estimate from the data presented in [ll[ that 
the laser intensity required to generate a Rabi frequency 
of 100 cm -1 is about 500 MW/cm 2 , which is comparable 



to the intensities used in standard CARS experiments, 
and sufficiently low to avoid multiphoton excitation [l2] ] . 
However, it should be noted that such a value will depend 
on the molecular system under study. 



IV. CONCLUSIONS 

In conclusion, we have demonstrated a route to obtain- 
ing sub-diffraction-limited CARS signals. By resonantly 
driving the excitation of a control state with relatively 
long decoherence times, intensity and thus spatially de- 
pendent sidebands in the CARS emission spectrum can 
be generated. By accurately measuring this Rabi split- 
ting in the CARS spectrum, objects can be resolved to a 
resolution below the diffraction-limit. Using typical num- 
bers, we show that a resolution in the order of 65 nm 
may be achievable. We note that our calculations also 
show the generation of Rabi sidebands in other molecule- 
specific emission lines, such as Coherent Stokes Raman 
Scattering. The approach towards sub-diffraction-limited 
resolution presented here can thus be used with this mi- 
croscopy technique as well, making it more broadly ap- 
plicable. 
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